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The Relativistic Heavy Ion Collider collides heavy nuclei at ultrarelativistic energies, creating a strongly interact- 
ing, partonic medium that is opaque to the passage of high energy quarks and gluons. Direct jet reconstruction 
applied to these collision systems provides a crucial constraint on the mechanism for in-medium parton energy 
loss and jet-medium interactions. However, traditional jet reconstruction algorithm operating in the large soft 
background at RHIC give rise to fake jets well above the intrinsic production rate of high-p T partons, imped- 
ing the detection of the low cross section jet signal at RHIC energies. We developed a new jet reconstruction 
algorithm that uses a Gaussian filter to locate and reconstruct the jet energy. This algorithm is combined with 
a fake jet rejection scheme that provides efficient jet reconstruction with acceptable fake rate in a background 
environment up to the central Au + Au collision at mn — 200 GeV. We present results of its application in 
p + p and Cu + Cu collisions using data from the PHENIX detector, namely p + p cross section, Cu + Cu jet 
yields, the Cu + Cu nuclear modification factor, and Cu + Cu jet-jet azimuthal correlation. 



1. Introduction 

Measurements of single particle production at the 
Relativistic Heavy Ion Collider (RHIC) have found 
significant suppression in Cu + Cu and Au + Au colli- 
sions at tJsnn — 130 GeV and y/s^N — 200 GeV (e.g. 
[l], Q). However, single or few-particle observables are 
fragmentation dependent, and only indirectly probe 
the energy loss of the parent parton. While direct 
jet reconstruction has become widely used to study 
perturbative quantum chromo dynamics (PQCD) in 
e + e" and hadronic colliders, the high multiplicity and 
strongly fluctuating background in heavy ion collisions 
made direct applications of jet reconstruction difficult. 

In PHENIX, the limited detector aperture presents 
additional difficulty to directly apply traditional jet 
reconstruction algorithms that are sensitive to large 
angle fragments (or the lack of such), such as the 
cone [1] and k± algorithms 0, [f| . PHENIX has two 
"central arm" spectrometers with an aperture consist- 
ing of two "arms" with \rj\ < 0.35 and Atfi = ir/2 
each Q. Since PHENIX is designed to measure rare 
and electromagnetic probes, and therefore can sus- 
tain a high read-out rate, it is well-suited to mea- 
sure high-p T jets at a low cross section. However, 
the narrow PHENIX central arm aperture especially 
in An leads to the loss of large angle fragments 
from jets and potentially larger systematic errors due 
to the importance of edge effects. With increasing 
Pt, the fragments of jets become increasingly colli- 
mated, with most of the energy concentrated in a 
cone much smaller than the PHENIX acceptance. A 
jet reconstruction algorithm that emphasizes the core 
and deemphasizes the large angle tail would therefore 
be less sensitive to both the background fluctuation 
in heavy ion events and the limited aperture of the 
PHENIX detector. 

In order to provide an effective method to recon- 
struct jets at the presence of heavy ion background 
and/or limited detector aperture, we started in 2006 



to develop a jet reconstruction algorithm that takes 
advantage of the collimated emission of hadrons and 
is insensitive to large angle fragments. We observed 
that the flat weighting in traditional jet reconstruction 
algorithms makes it prone to large angle fluctuations. 
While the background grows with R 2 , the jet contri- 
bution to the pt grows slowly above R > 0.3 (e.g. Q)- 
A nonflat weighting that enhance the core jet signal 
to periphery would naturally suppress this sensitivit 
and stabilize the jet axis. The energy flow variable 
also provided us with the hint that angular convolu- 
tion of the event transverse momentum (p T ) density 
with a continuous distribution can be an effective de- 
scription of QCD processes. We are therefore using a 
Gaussian filter as a generalized form of the cone algo- 
rithm to reconstruct jets (d, 03). Early in the devel- 
opment of jet reconstruction algorithms, the British- 
French-Scandinavian collaboration used the equiva- 
lent to a Gaussian filter with o = 0.5 [ll|, preceding 
even the Snowmass accord on the cone algorithm Q . 

The heavy ion background itself is generated by 
combination of soft, collective processes, semihard 
processes and subsequent QCD hadronization. The 
presence of heavy flavor decays and minijet produc- 
tion suggest that one should expect the background to 
contain more complex, angularly correlated structures 
than from a purely randomized and isotropic process. 
Therefore, a strategy based on statistically subtract- 
ing the background would require detailed knowledge 
about the structure of the heavy ion background and 
its interaction with the jet reconstruction algorithm. 

Fluctuations in the underlying event of heavy ion 
collisions arc known to cause a false apparent jet pro- 
duction, if a jet reconstruction algorithm is applied 
without proper handling of the background 12]. At 
the collision energy of RHIC, the intrinsic jet pro- 
duction rate is far below the high-p T tail of the 
background fluctuation across a large p T range that 
is statistically accessible to RHIC. Based on the 
p + p jet cross-sections presented later in this pa- 
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per, we estimate that for the central Cu + Cu at 
^/sWn = 200 GeV and p T between 10-20 GeV/c and 
20-30 GeV/c, the jet yields are N^dN/dy w 4 x 10" 3 
and 8 x 10~ 5 , respectively, which means that rare fluc- 
tuations can contribute significantly to the jet yield, 
if the background cannot be suppressed below this 
level. Avoiding the px range with background contri- 
bution would require measuring jets only at the very 
end of the pr range, where the statistics are poor. 
The presence of rare fluctuations also suggest that a 
simple removal of \ow-p T particles would only sam- 
ple the stronger fluctuations in the rare, high-p T tail, 
while biasing both the fragmentation and energy of 
the reconstructed jets. 

We therefore developed a fake jet rejection strategy 
that is based on the jet versus background shape and 
can achieve a higher rejection rate than previously 
proposed algorithms for the LHC (e.g. 13]). As will 
be demonstrated below, the fake rejection provides a 
fast rise to unity efficiency within the RHIC accessible 
p T range. 



2. Jet reconstruction by Gaussian 
filtering 

It can be shown that the iterative cone algorithm is 
equivalent to finding local maxima of a filter output in 
(?7, <f) with a flat angular weighting k(r 2 ) = 0(R 2 — r 2 ) 
with r 2 = rj 2 + <j) 2 (note that unlike k{r 2 ), the filter 
kernel h(r 2 ) oc — J dr 2 k{r 2 ) oc max(0, 1 — r 2 /R 2 ), and 
not flat) [HOI. The cone algorithm entails a spe- 
cific choice of angular weighting. The Gaussian filter 



is based on another, which takes advantage of jets be- 
ing a collimated emission of particles, and enhances 
the center of the jet and suppresses the possible con- 
tribution from the event background in the periphery. 
Expressed in the filter form described below, the al- 
gorithm samples the entire possible (rj, <ft) range and 
is seedless. By additionally avoiding a sharp radial 
cutoff, the algorithm therefore becomes analytically 
collinear and infrared safe (we also verified the practi- 
cal infrared safety using a procedure similar to [HI). 

A combined event transverse momentum density 
that contains both the final state particles pT,i, and 
p-]?(T],(f>)) which represents an independently evalu- 
ated average contribution from the underlying event, 
can be defined as 

Pt(??» = ^2pt,iS(v ~ m)$(<t> - </>i) - P^iihfi), (!) 

i(EF 

The Gaussian filtering of p T is the linear-circular con- 
volution of pt(Vi <$) with a Gaussian distribution 

PtXv,4>) =JJ ^'#V(77',0')e- ((r '- r '' )2+( ^ ' )2)/2CT . 

The output of the filter for a given (rj, <f) position is 
the Gaussian-weighted transverse momentum in that 
event above the average background from the under- 
lying event. The local maxima in p T lt (r], </>) are the 
reconstructed jets using the Gaussian filter. 

Figures HH2] demonstrate the behavior for a pip 
and Cu + Cu event, respectively. Charged tracks and 
photons are shown at the bottom by a Lego plot. 
The distribution of filter output values of the event 
is shown at the top as a contour plot. The maxima in 




Figure 1: A PHENIX Run- 5 p + p at yfs = 200 GeV di- 
jet event. Charged tracks and photons are shown at the 
bottom by a Lego plot. The distribution of filter output 
values of the event is shown at the top as a contour plot. 
The maxima in the filter density are reconstructed as jet 
axes, shown as red lines at the positions on the contour 
and Lego plots. 




Figure 2: A PHENIX Run-5 Cu + Cu at ^i^" = 200 GeV 
dijet event at ~ 20% centrality. Charged tracks and pho- 
tons are shown at the bottom by a Lego plot. The distri- 
bution of filter output values of the event is shown at the 
top as a contour plot. The maxima in the filter density 
are reconstructed as jet axes, shown as red lines at the 
positions on the contour and Lego plots. 
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Figure 3: Comparison of the jet energy scale between the 
a = 0.3 Gaussian filter (p^ H ) and the R = 0.4 SISCone 
(p| ISCone , with 0.5 overlap threshold) [lifer p + p col- 
lisions at y/s — 200 GeV using PYTHlAyjJ and central 
Au + Au collisions at ^s NN = 200 GeV using pyquench 
[lc| at the event generator level (no detector effects). The 
dashed line indicates the position of pf. 1 * 
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the filter density are reconstructed as jet axes, shown 
as red lines at the positions on the contour and Lego 
plots. 

The angular weighting of the Gaussian filter modi- 
fies the energy summation and therefore will not pro- 
duce the same jet energy as e.g. the cone algorithm. 
To evaluate the difference between the energy scales of 
the a = 0.3 Gaussian filter and the R = 0.4 SISCone 
algorithms (0.5 overlap threshold) [lq . we show in 
Figure[3]the jet energy scale comparison for both p + p 
collisions at yfs — 200 GeV using pythia [ijj and 
central Au + Au collisions at ^Jsmn = 200 GeV using 
PYQUENCH [l8| at the event generator level. Further 
study of its property in p + p events are detailed in 

In heavy ion events, the background contribution 
p T s will depend on the collision centrality and reac- 
tion plane. The Gaussian filter allows these collision 
variables to be fully parameterized, which is crucial 
for heavy ion jet reconstruction algorithm. Due to 
the narrow pseudorapidity coverage at PHENIX, we 
also parametrize depending on the collision vertex 
z position. Since in the Cu + Cu system at ^/snn = 
200 GeV, the event-by-event fluctuation strongly dom- 
inates over flow effects, we do not parametrize with 
respect to the reaction plane (note that the linearity 
in the jet definition means that a weak flow simply 
translates into an additional p T smearing). 




Figure 4: The PHENIX central arm detectors for RHIC 
Run-5 (year 2004/2005), viewed along the beam axis from 
the south towards north. Dark regions indicate detectors 
used for the jet reconstruction: The drift chamber (DC), 
the pad chamber layers 1 and 3 (PC1/PC3), the ring- 
imaging Cerenkov detector (RICH) , and the Pb scintillator 
(PbSc) and Pb glass (PbGl) electromagnetic calorimeters. 



central arm detectors used for jet reconstruction are 
the drift chamber (DC), the pad chamber layers 1 
and 3 (PC1/PC3), the ring-imaging Cerenkov de- 
tector (RICH), and the electromagnetic calorimeters 
(EMCal). For the data presented in this paper, 
DC/PC1/PC3 provide momentum measurement for 
charged particles, and the EMCal the energy for pho- 
tons and electrons. Two calorimeter technologies 
were used, 6 of the total 8 sectors are covered by 
Pb-scintillator (PbSc) EMCal, 2 sectors by Pb-glass 
(PbGl) calorimeters. 

Pattern recognition and momentum reconstruction 
of the tracking system is performed using a combina- 
torial Hough transform. The px scale is determined 
by time-of-flight measurement of ^ , K , and p/p. 

Since PHENIX currently does not have the ability 
of performing in-field tracking 1 , conversion electrons 
in the DC can produce a displaced track that appears 
to the momentum reconstruction to be a very high- 
ly track originating from the event vertex. Informa- 
tion from the RICH and the dE/dx measurement is 
therefore used to identify and remove these conversion 
electrons. To provide additional suppression for jets 
with > 20GeV/c, we use the fact that conver- 
sion electrons are unlikely to angularly coincide with 
a high-p T jet, and require the reconstructed jet to 
have a minimum jet multiplicity of 3 particles, and 
the charged fraction of the jet p T to be below 0.9 to 
remove events with single track and falsely large px 
values. 



3. Experimental setup 

FigureUshows the Run-5 I uumiv ^uuoi aim ^uu- cap 
figuration for RHIC Run-5 (year 2004/2005). The [3 



x a Si vertex detector upgrade is going to provide such a 
FigureHshowS the Run-5 PHENIX central arm COn- capability and is scheduled to be installed in year 2010, see e.g. 
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Shower shape cuts were applied to the EMCal clus- 
ters to remove clusters generated by hadronic showers. 
The absolute energy scale of the calorimeter clusters 
is determined both by the reconstructed ir° masses 
from 7T° — » 77 decay, and checked by comparing RICH 
identified momenta from tracking against the cor- 
responding cluster energies. The residual uncertainty 
in the energy scale is ±3% (syst.). 

The PHENIX minimum bias (MB) trigger is defined 
by the coincident firing of the two beam-beam coun- 
ters (BBC) located at 3.0 < r\ < 3.9. The Run-5 BBC 
cross section is ctbbc = 22.9 ± 2.3 mb (syst.), mea- 
sured using the Van de Meer/ vernier scan method. 
The efficiency of BBC firing on an event containing a 
jet with p^ c > 2 GeV/c is £ B bc = 0.86 ± 0.05 (syst.) 
and within that uncertainty, approximately constant 
with respect to p™ c . For both our p + p and Cu + Cu 
measurements, we require the collision vertex to be 
within \z\ < 25cm along the beam axis, derived from 
the BBC timing information. 

All the charge tracks and electromagnetic clusters 
in the EMCal passing the described cuts are used in 
the Gaussian filter, and all resulting maxima are con- 
sidered candidate jets. Since we do not explicitly split 
jet pairs with small angular separation, it is possible 
to reconstruct jet pairs with substantial overlap. In 
reality however, we rarely observe jets reconstructed 
with an angular separation of AR < 0.5. 



4. Jet spectrum in p + p 

The data presented in this section were obtained 
from the PHENIX p + p dataset from the RHIC Run- 
5 (year 2004/2005). After removal of bad quality runs, 
a total of 1.47 x 10 9 minimum bias p + p and 1.16 x 10 9 
triggered p + p events are being used. 

PHENIX can trigger on high-p T and electromag- 
netic processes using the central arm EMCal-RICH- 
trigger (ERT). For the result presented here, we use a 
trigger that requires a total energy E > 1.4 GeV de- 
posited in a 4 x 4 group of calorimeter towers. This 
trigger is well-suited for jet measurement due to its 
low noise level and fast efficiency saturation with re- 
spect to the jet p T . The efficiency is evaluated in term 
of the p 1 ]; c of the most energetic jet in the event, and 
rises from approximately 0.25 for p^? c 2.2 GeV/c 
to 0.9 for p^ c ~ 9.5 GeV/c. After correction for the 
efficiency, the minimum bias and triggered datasets 
agree within an uncertainty of 5% (syst.) in the range 
from p 1 ^ = 2GeV/c up to 30GeV/c (limited by the 
minimum bias data set statistics). Therefore, we cal- 
culate a combined jet yield per minimum bias event 
according to 

1 dN _ 1 £ ertU>t ) dp»» + Ip^T 
iVevt dp 1 ? N M B smb + 1 




Figure 5: The PHENIX jet P{p^\p T ) transfer matrix for 
V^s = 200 GeV and a — 0.3 Gaussian filter, derived from 
the GEANT simulation of w 1.6 x 10 7 pythia events. The 
V V t < Vt region is dominated by n, K\ energy loss. 

where smb = 37.34 is the average scale down applied 
to the minimum bias trigger. 

The spectrum is then unfolded by the regularized 
inversion of the Fredholm equation 

dN f , / w i \ dN 
^ = J dPTP(f T *\ PT )^ (4) 

using singular value decomposition (SVD). The si- 
multaneous minimization of the second order finite- 
differences is use as a constraint to the continuity of 
the unfolding result. This procedure is implemented 
in the software package GURU [2(|. The transfer ma- 
trix of P{p p ^\pt) expresses the (conditional) proba- 
bility that a jet with the true transverse momentum, 
p Tl is reconstructed with p^ . 

We used pythia 6.4.20 with the (6.4.20-default) 
"old" multiparton interaction scheme and GEANT sim- 
ulation to evaluated P(p^ \pt)- A total of 1.6 x 10 7 
events were simulated with 1 4 d ifferent minimum 
Q 2 settings varying betweeni/Q^ > 0.5GeV/c and 
y/Q 2 ' > 64GeV/c. The transfer matrix Pij^lpr) 
resulting from the simulation procedure is shown in 
figure [5] The p^ < pt region is dominated by n, 
energy loss. 

We measured the y/s = 200 GeV p + p spectrum 
using the combined PHENIX Run-5 minimum bias 
and triggered data. We require that the reconstructed 
jet is at least by an angular distance of d > 0.05 rad 
inside the infinite momentum PHENIX central arm 
acceptance. The spectrum with respect to the true 
Pt is then given by 

Ed 3 a 1 dPu _ ctbbc I I dN 
dp 3 2ttpt dprdy A 6bbc Pt N cvt dpr 

where dN/ dpx is the unfolding result from and 

A= 2(Ar]-2d)(A(j)/2-2d) (6) 
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Figure 6: PHENIX Run-5 p + p at y/s = 200 GeV invari- 
ant jet cross section spectrum as a function of pr- The 
shaded box to the left indicates the overall normaliza- 
tion systematic uncertainty, shaded boxes associated with 
data points indicate point-to-point systematic uncertain- 
ties, and error bars indicate statistical uncertainties. 

the fiducially reduced PHENIX central arm accep- 
tance area. 

The regularized least square unfolding involves the 
regularization parameter r (or sometimes A 2 ) and its 
choice translates into an uncertainty on the global 
shape/low frequency component of a jet spectrum. 
We evaluated the systematic uncertainty in the un- 
folded spectrum due to such variation in the regular- 
ization parameter by varying r over the entire mean- 
ingful range between « 4 degrees of freedom up to the 
Nyquist frequency. We combine the resulting, point- 
by-point estimate of the systematic uncertainty for 
the spectrum shape as part of the total experimental 
systematic uncertainty. The so evaluated systematic 
uncertainty is therefore representative of the full range 
of regularization parameter choices. 

The residual ±3% systematic uncertainty in the en- 
ergy scale translates into a constant uncertainty of 
±15% (syst.) for the jet spectrum (due to its power- 
law shape with largely constant exponent). 

Figure [5] shows the PHENIX preliminary p + p jet 
spectrum measured using the Gaussian filter by the 
procedures described above, plotted in invariant cross 
sections. The shaded box to the left indicates the 
overall normalization systematic uncertainty, shaded 
boxes associated with data points indicate point-to- 
point systematic uncertainties, and error bars indicate 
statistical uncertainties. We show the unfolded spec- 
trum out to the pt bin where the nominal yield for 
the number of sampled events reaches the level of 1 
jet, namely 60GeV/c. 

Figure [7] shows the same spectrum as in Figure [6l 
compared against the spectrum from [2l| . the next- 
to-leading order (NLO) calculation using the small 
cone approximation (SCA) (22} , and the leading order 



Figure 7: PHENIX Run-5 p + p at y/s = 200 GeV invari- 
ant jet cross section spectrum as a function of pr, with 
comparison to [21|, next-to- leading order calculation from 
[22l ]. and pythia assuming K = 2.5. The shaded box 
to the left indicates the overall normalization systematic 
uncertainty, shaded boxes associated with data points in- 
dicate point-to-point systematic uncertainties, and error 
bars indicate statistical uncertainties. 

PYTHIA spectrum assuming K = 2.5. The comparison 
to [2l| and NLO SCA involve different jet definitions, 
a residual difference should be expected, even though 
for px > 15 GeV/c it appears to be small between fil- 
ter and cone jets for the Gaussian size rj = 0.3 used in 
this analysis. Our spectrum is close to [2l[ within its 
p T reach. The spectrum also follows approximately 
the shape of the NLO SCA calculation, and the lead- 
ing order pythia spectrum, if K = 2.5 is assumed. 
However, a more appropriate comparison would in- 
volve Gaussian filter based NLO calculations, which 
we plan to perform in the future. 



5. Cu + Cu results 

The data presented in this section were obtained 
from the PHENIX Cu + Cu dataset from the RHIC 
Run-5 (year 2004/2005). The p + p data presented 
above provide a baseline for the Cu + Cu measure- 
ments and explicitly appear in the jet nuclear modifi- 
cation factor. After removal of bad quality runs, a to- 
tal of 1.58 x 10 8 minimum bias are used in the Cu + Cu 
analysis, covering the centrality range 0-95%. Be- 
cause of poor statistics and large uncertainties in the 
Tab for very peripheral events, we excluded events in 
the 80-95% centrality range and divide the remainder 
into 4 bins. 

The p T s in |T]) is parametrized as a product of the 
centrality and vertex dependent total event t times 
a vertex dependent priv^) /Pt* distribution. The 
values are estimated by averaging over the minimum 
bias Cu + Cu events. 
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Figure 8: The a = 0.3 Gaussian filter p T resolution with 
Apr /pt oc 1 /pr (constant energy smearing) fits for differ- 
ent centralities of Cu + Cu collisions at ^/snn = 200 GeV, 
evaluated by embedding PHENIX Run-5 p + p events into 
PHENIX Run-5 Cu + Cu events. The constant energy 
smearing behavior is expected for an (asymptotically) lin- 
ear jet reconstruction algorithm. 



In heavy ion collisions, centrality-dependent fluctu- 
ations in the underlying event will broaden the energy 
resolution of the jet measurement beyond the detector 
response described in Section 4. The differences be- 
tween the Cu + Cu jet p T scale, p^ uCu and the p + p 
energy scale must be accounted for in any comparison 
of the jet spectra or evaluation of Raa- 

We evaluate the effects of the Cu + Cu underlying 
event on jet measurements by embedding PHENIX- 
measurcd p + p events into the minimum bias Cu + Cu 
events by combining their list of tracks and clusters. 
We reconstruct the embedded events and evaluate 
the energy and angular resolution and the reconstruc- 
tion efficiency for jets that match the input p + p jets 
within a radial separation AR < 0.3. Both events 
are required to fall into the same Az = 5 cm bin 
in vertex position. The evaluated transfer matrix is 
then recombined into the final centrality binning us- 
ing AT co n scaling. Due to the linearity of Gaussian 
filter, for p^f > 8GeV/c the energy smearing is con- 
stant within our statistical uncertainty. We use this 
feature of the energy response to extend the trans- 
fer matrix to higher p T where the limited statistics 
in the p + p sample would otherwise introduce fluc- 
tuations into the transfer matrix that would have a 
negative impact on the unfolding. 

Figure [8] shows the a = 0.3 Gaussian filter p T 
resolution for different centralities of Cu + Cu colli- 
200 GeV, evaluated by embedding 
■p events into PHENIX Run-5 



sions at ^/snn 
PHENIX Run-5 p 
Cu + Cu events. Superimposed are Apx/pr oc l/pr 
fits that shows the constant energy smearing, which is 
expected for an (asymptotically) linear jet reconstruc- 
tion algorithm such as the Gaussian filter. 




Figure 9: Comparison between a PHENIX Run-5 Cu + Cu 
event with a 9.6GeV/c jet passing fake rejection and an 
event with 10.8GeV/c "jet" that fails the fake rejection, 
and is presumed to be a background fluctuation 



Due to the large N co \\ in central Cu+Cu events, 
the underlying event has potentially a high yield to 
contain an intrinsic jet. This is avoided by requir- 
ing that reconstructed jet after embedding matches a 
p + p jet with p^P > 4 GeV /c within an angular range 
of AR < 0.3. The residual contamination was evalu- 
ated by reconstructing the underlying event jet spec- 
trum within AR of the original p + p jet axis, and is 
found to be < 10~ 2 . 

In the measurement presented here, we use two 
approaches to correct for the energy scale difference 
in Cu+Cu collisions: (a) by unfolding the Cu + Cu 
spectrum using the transfer matrix P{pt uCu \Pt')j the 
Raa is then calculated by comparing against the p + p 
spectrum; and (b) by embedding p + p events into 
Cu + Cu events, so that the p + p spectra attains the 
Cu + Cu background induced energy smearing, and 
the Raa is derived by comparing the Cu + Cu spec- 
trum to the embedded spectrum. 



5.1 . Fake jet rejection 

In order to reduce the contribution of fake jets, we 
apply a shape based discriminant to the jets. The 
discriminant is defined as 



9<T dl s (V, </>) 



E 

i£ fragment 



Pr,i e 



-v) 2 +(<l>i 



(7) 

where (77, <fi) is the reconstructed jet axis. The dis- 
criminant size Odis = 0.1 to chosen to be approx- 
imately or below the size AR h s = y/2ir/{dN/dT]), 
which is the characteristic background particle sep- 
aration. If the local distribution of high-p T parti- 
cle indicates that the jet may be misreconstructed in 
its angle (e.g. a background fluctuation beneath the 
jet with a nonzero gradient), we also search for an 
(77, 4>) where <7<r dia might be maximized. We denote 
this "adapted" discriminant g' . In hijing studies, 
we found this discriminant, when applied to central 
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Figure 10: Efficiency for the g' 01 > 52(GeV/c) 2 fake jet 
rejection that is suitable for central Au + Au collisions at 
^snn ~ 200 GeV for p + p collisions at y/s — 200 GeV us- 
ing pythia [ItJ and central Au + Au collisions at ^/snn = 
200 GeV using pyquench [l]| at the event generator level 
(no detector effects), compared to the pythia efficiency for 
g'o.i > 17.8 (GeV/c) 2 that is used for Cu + Cu collisions. 



Au + Au collisions at RHIC, to significantly outper- 
form discriminants proposed for the LHC, such as the 
Pt/{A) [13 or Ej T 0. 

Figure[9]illustrates the behavior of the fake jet rejec- 
tion using two actual PHENIX Run-5 Cu + Cu events. 

With the definition in ((?]), it is also possible to 
achieve a largely centrality independent efficiency 
turn-on. This feature is important in order to con- 
strain the systematic impact in a centrality dependent 
measurement. We could, in principle, adapt the fake 
rejection to the collision centrality but that by keeping 
the rejection threshold fixed we obtain a nearly cen- 
trality independent jet finding efficiency. Using the 
studies of the discriminant distribution in Cu + Cu 
collisions and on the dijet A(j) analysis presented be- 
low, we have chosen a nominal discriminant threshold 
oig' 01 > 17.8(GeV/c) 2 . 

There are several methods to investigate the sen- 
sitivity of reconstructed jets to the fake rejection. 
The discriminant distribution shows that for px > 
16GeV/c, the jets have discriminant values that lie 
mostly above our Cu + Cu fake jet rejection threshold 
of g' > 17.8 (GeV/c) 2 . Changing or removing the 
threshold therefore has little impact on our spectra 
and Raa above 12-16 GeV/c. 

We further studied the sensitivity of both efficiency 
after fake jet rejection and the energy scale using 
PYTHIA and PYQUENCH at the fake jet rejection level 
<7o j > 52 (GeV/c) 2 required for central Au + Au col- 
lisions at ^/sjvn — 200 GeV. Both the efficiency 
and the jet energy scale compared to SISCone shows 
very small modification, i.e. below the present sys- 
tematic uncertainty in the overall energy scale. This 
is shown in Figure 1101 and the comparison between 
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Figure 11: Efficiency for the g' 0A > 17.8 (GeV/c) 2 fake 
jet rejection for different centralities of Cu + Cu collisions 
at v /i^7 = 200 GeV, evaluated by embedding PHENIX 
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Figure 12: Run-5 Cu + Cu at -y/ijvjv = 200 GeV azimuthal 
correlation for 0-20% centrality dijets in yields and with 
different g'o.i fake jet rejection thresholds. The nominal 
fake jet rejection threshold used is g' i > 17.8 (GeV/c) 2 . 



pythia and pyquench jet energy scales was shown 
in Figure [31 Figure QT] shows the efficiency for the 
g'o i > 17.8 (GeV/c) 2 fake rejection for different cen- 
tralities of Cu + Cu collisions at ^snn = 200 GeV. 

Angular correlation of dijets can be used to study 
the yield of residual fake jets. The yield of three and 
more jets at large angle with respect the leading dijet 
axis (e.g. due to multiple collisions or initial and final 
state radiation) is strongly suppressed at high-pr, and 
therefore the decorrelated yield approximately per- 
pendicular to the dijet axis is a good estimator of the 
residual contamination by fake jets. In the Cu + Cu 
data, the A</> distribution saturates at our normal 
discriminant threshold of g' Q x > 17.8 (GeV/c) 2 . At 
p£uCu _ 7.5GeV/c for 0-20% centrality, this trans- 
lates into an upper bound of 10% for the fake jet con- 
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Figure 13: PHENIX Run-5 Cu + Cu at ^sWn = 200 GeV 
invariant jet yield as function of p^P, with comparison to 
the {Tab} scaled p + p cross section. The shaded box to 
the left indicates the centrality-dependent systematic un- 
certainty in the normalization, shaded boxes associated 
with data points indicate point-to-point systematic uncer- 
tainties, and error bars indicate statistical uncertainties. 



Figure 14: PHENIX Run-5 Cu + Cu at ^/sITn = 200 GeV 
Raa derived from unfolding. The shaded box to the left 
indicates the p + p-Cu + Cu systematic uncertainty in the 
jet energy scale, shaded boxes to the right shows central- 
ity dependent systematic uncertainty between embedding 
and unfolding, shaded boxes associated with data points 
indicate point-to-point systematic uncertainties, and error 
bars indicate statistical uncertainties. 



tamination in the jet yield. Figure [T2l shows the effect 
of the increasing fake jet rejection on the A(f> distribu- 
tion using the 0-20% centrality Run-5 Cu + Cu data. 

5.2. Cu + Cu jet spectra and nuclear 
modification factors 

We unfold the raw spectrum dN / dp? using the 
centrality dependent transfer matrix P(px uCu \Pt') 
and obtain the spectrum in term of the dp^F energy 
scale, dN/dp^?. The invariant yield in the respective 
Pt scale is then given by 

1 Ed 3 N _ 1 1 1 1 dN 
TVcvt dp 3 e b (p T ) Ap T iVe Vt dpT 

with the same fiducially reduced area A defined in ^ , 
and where £f r is the fake rejection efficiency (compare 
figure [TTj). For both the jet spectra and Raa measure- 
ment, we conservatively restrict our pt range to the 
region with ef r > 0.75. When applying the charged 
fraction cut in Cu + Cu events, we only include parti- 
cles with pt > 1.5 GeV/c, so that background fluctu- 
ations cannot cause us to accept fake jets generated 
by late conversion electrons. 

Figure [T2] shows the invariant jet yield as function 
of p V r^ together with the (Tab) scaled p + p cross sec- 
tion, where Tab is the nuclear overlap function with 
Cu + Cu values given in are used. The shaded box 
to the left indicates the centrality-dependent system- 
atic uncertainty in the normalization, shaded boxes 
associated with data points indicate point-to-point 
systematic uncertainties, and error bars indicate sta- 
tistical uncertainties. 



The nuclear modification factor is defined by 

Raa = N <^ dNcnCn / dpT (9) 
(T A B)da pp /dp T 

For the Raa in V V t ener gy scale and using unfolding, 
we divide the spectra shown in Figure [T31 by the raw 
p+p spectrum in p^P . 

To check for potential systematic errors resulting 
from the unfolding, we alternatively evaluated the 
Raa by comparing the raw Cu + Cu spectrum to the 
embedded p + p spectrum. Since the extracted Raa 
do not show a significant p T dependence within the 
uncertainties, we can compare between the two sets 
of Raa despite the difference in the jet energy scale. 
We additionally studied the sensitivity of the Raa to 
vertex and fiducial cuts. The so evaluated, centrality 
dependent differences from the unfolding/embedding 
comparison, vertex dependence, fiducial dependence, 
are combined to obtain the centrality dependent sys- 
tematic uncertainties. 

Figure Q3] shows the extracted Raa using unfold- 
ing. The shaded box to the left indicates the p + p- 
Cu+Cu systematic uncertainty in the jet energy 
scale, shaded boxes to the right shows centrality de- 
pendent systematic uncertainty between embedding 
and unfolding, shaded boxes associated with data 
points indicate point-to-point systematic uncertain- 
ties, and error bars indicate statistical uncertainties. 

Figure [T5] shows the comparison between the Raa 
derived using unfolding and embedding (with the 
same notation as in Figure [P4")) . The level of the 
suppression obtained with both methods is the same 
within the statistical and unfolding-systematic uncer- 
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Figure 15: Comparison between the PHENIX Run-5 Cu + Cu at ^/sNN = 200 GeV Raa derived from unfolding (filled 
symbols) and embedding (open symbols). The shaded box to the left indicates the p + p-Cu + Cu systematic uncertainty 
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Figure 16: Comparison between the central PHENIX 
Run-5 Cu + Cu at v / sjvjv = 200 GeV jet Raa derived from 
unfolding and the 7r° Raa ■ The shaded box to the left in- 
dicates the p + p-Cu + Cu systematic uncertainty in the 
jet energy scale, shaded boxes to the right shows central- 
ity dependent systematic uncertainty between embedding 
and unfolding, shaded boxes associated with data points 
indicate point-to-point systematic uncertainties, and er- 
ror bars indicate statistical uncertainties. Note that while 
the flatness of Raa makes a comparison across different 
energy scales possible, 7r° with (z) = 0.7 has a different 
energy scale. 



Centrality 


Width 


0-20% 


0.223 ±0.017 


20-40% 


0.231 ± 0.016 


40-60% 


0.260 ± 0.059 


60-80% 


0.253 ±0.055 



Table I Widths of Gaussian fit to the PHENIX Run-5 
Cu ± Cu at y / sjvjv = 200 GeV azimuthal angular correla- 
tion for jets with 7.5 GeV/c < p£ uCu < 11.5 GeV/c 



therefore allows a comparison. 

We observe a Raa that becomes gradually sup- 
pressed with increasing centrality. The level of sup- 
pression in the most central 20% centralities is at 
Raa ~ 0.5-0.6 and comparable to that of ir°. 



5.3. Cu + Cu jet-jet azimuthal correlations 

The Cu + Cu jet-jet azimuthal correlation is ex- 
tracted by correcting for the acceptance effect using 
the area-normalized mixed event yield (e.g. [HI): 



dN{A(j)) _ 1 dN Iaw (A4>) 



(10) 



tainties, which gives us confidence that the unfolding 
procedure is not significantly biasing the result. Fig- 
ure [TH] compares the central 20% suppression with the 
7T° suppression from [2j (with the same notation as in 
Figure [T5f . While the Raa of tt° has a different en- 
ergy scale than jets, both Raa are approximately flat 
with respect to p T within our accessible range and 



where A(A<f>) is the detector acceptance correction. 
Using a Gaussian fit to the distribution, we extracted 
the width for 7.5GeV/c < p£ uCu < 11.5GeV/c. The 
widths are consistent within the uncertainty across all 
centrality ranges. 

Figure [17] shows the azimuthal jet-jet correlation 
with Gaussian fits for jets with 7.5 GeV/c < p£ uCu < 
11.5GeV/c. Table [T] lists the Gaussian widths ex- 
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Figure 17: PHENIX Run-5 Cu + Cu at </sWn = 200 GeV 
azimuthal jet-jet correlation with Gaussian fits for jets 
with 7.5 GeV/c < p? uCu < 11.5 GeV/c 



tracted from the azimuthal jet-jet correlation for jets 
with 7.5 GeV/c < p£ uCu < 11.5 GeV/c. 



6. Discussion 

Using the PHENIX Run-5 p + p dataset, we ex- 
tracted the first RHIC p+p spectrum that extends to 
p T « 60GeV/c or Ed 3 a/dp 3 w 0.3 fb GeV~ 2 c 3 . This 
effectively demonstrates the capability of PHENIX as 
a detector for the study of jet physics. 

Applying our algorithm to Cu + Cu collisions, were 
also able to clearly demonstrate, for the first time, 
the feasibility of direct jet reconstruction in heavy ion 
collisions and for all centralities. We are able to show 
the centrality dependent onset of jet suppression from 
peripheral to central collision. This demonstrates the 
capability of the Gaussian filter algorithm as a heavy 
ion jet reconstruction algorithm, that can be applied, 
fully unconstrained, in any heavy ion collision setting 
at RHIC. 

Our a = 0.3 Raa in Cu + Cu suggest that the jet 
production is strongly suppressed, and is compara- 
ble to the suppression level of high-p T 7r°. Strong 
suppression in reconstructed jets indicate that signif- 
icant amount of energy disappears from the angular 
region covered by the size of the jet reconstruction 
algorithm. Since the fake jet rejection has little im- 
pact on jets with pt > 12-16 GeV/c, we observe an 
overall suppression that is independent from the fake 
rejection scheme. However, if jets significantly broad- 
ens at the presence of the medium, it is possible that 
small angle jet reconstruction algorithm such as the 
a = 0.3 Gaussian filter is deselecting these jets. We 
are current pursuing jet spectra and Raa evaluation 
with larger angle to address this. 

The azimuthal correlation suggest that the sur- 
viving parton traversing the medium has very small 



transverse kx broadening. 

While our current results are encouraging and 
demonstrates the accessibility of a large inventory of 
heavy ion jet variables to PHENIX, we would have to 
continue our measurement in term of fragmentation 
properties and jet angular size dependence to under- 
stand and complete our picture of the parton energy 
loss. Since our p + p results probes in a previously 
unreached RHIC kinematic range, additional compar- 
isons with PQCD calculations could also provide valu- 
able insights. 
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